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ABSTRACT: The dumbbell structure of two-dimensional group IV material offers alternatives to grow thin films for diverse 
applications. Thermal properties are important for these applications. We obtain the lattice thermal conductivity of low-buckled 
(LB) and dumbbell (DB) silicene by using first-principles calculations and the Boltzmann transport equation for phonons. For 
LB silicene, the calculated lattice thermal conductivity with naturally occurring isotope concentrations is 27.72 W/mK. For DB 
silicene, the calculated value is 2.86 W/mK. The thermal conductivity for DB silicene is much lower than LB silicene due to 
stronger phonon scattering. Our results will induce further theoretical and experimental investigations on the thermoelectric 
(TE) properties of DB silicene. The size-dependent thermal conductivity in both LB and DB silicene is investigated as well for 
designing TE devices. This work sheds light on the manipulation of phonon transport in two-dimensional group IV materials by 


dumbbell structure formed from the addition of adatoms. 
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1. INTRODUCTION 


Unique electronic properties, high chemical stability, and 
mechanical strength have made graphene a material of interest 
in various fields ranging from electronics to energy conversion 
devices.” Reasonably, other group IV elements have been 
considered as candidates to have a graphenelike two-dimen- 
sional (2D) honeycomb structure with similar exceptional 
properties.” '* Recently, silicene, germanene, and stanene (2D 
silicon, germanium, and tin nanosheet, respectively) have been 
fabricated by epitaxial growth on substrates.'”-~’ Compared to 
other 2D group IV materials, the realization of silicene would 
be of great importance because of its easier incorporation into 
silicon-based microelectronics industry. First-principles calcu- 
lations have predicted that silicene has a low-buckled (LB) 
honeycomb structures.”'* However, to explain some character- 
istics of silicene grown on Ag(111) surfaces, some modifica- 
tions of this graphene-like structure have been proposed: the 
addition of Si adatoms to LB silicene results in the formation of 


a dumbbell (DB) structure with a lower total energy per 
atom.'°**"~* These works enrich the family of 2D group IV 
materials beyond the LB honeycomb lattices. Furthermore, DB 
structures occur not only for monolayer monatomic honey- 
comb structures, but also for their compounds such as SiGe and 
SiC. For the formation of DB structures, the energy 
barrier is not high, which can offer alternatives to grow thin 
films with diverse properties.” Although the electronic, optical, 
and magnetic properties of DB silicene have been studied 
intensively, %4 the phonon transport properties, which play 
an important role in potential applications such as electronic 
and energy conversion devices, are not revealed. 

For application in thermoelectric (TE) energy conversion, it 
is important to reduce the lattice thermal conductivity of a 


material while maintaining a high electrical conductivity, which 
often conflict with each other.’”** Traditionally there are four 
factors to guide the search for low thermal conductivity 
materials: (1) weak interatomic bonding, (2) high average 
atomic mass, (3) complex crystal structure, and (4) high 
anharmonicity.” According to Slack’s theory, weak interatomic 
bonding and high average atomic mass imply a low Debye 
temperature, and high anharmonicity indicates strong anhar- 
monic phonon scattering, which lead to a low thermal 
conductivity. This has been observed in Cu,GeSe, and 
monolayer transition metal dichalcogenides (MoS,, MoSe, 
and WS,) in our previous work.’**° Fully microscopic 
computational materials techniques have been developed to 
study thermal transport properties recently.” 7>? Thus, system- 
atic investigations of phonon transport properties in LB and 
DB silicene can be performed. 

Here we investigate the lattice thermal conductivity x of both 
LB and DB silicene using first-principles calculations and self- 
consistent iterative solution of the Boltzmann transport 
equation (BTE) for phonons.” *? A much lower thermal 
conductivity in DB silicene is observed (2.86 W/mK), which 
may indicate higher TE efficiency and will induce further 
investigations on the TE performance of DB silicene. The 
underlying mechanism is investigated by both Slack’s theory 
and calculation of phonon relaxation time. The influence of 
naturally occurring isotopes is studied as well. Finally, for 
designing TE nanostructures, we investigate the size-dependent 
lattice thermal conductivity. 


2. METHODOLOGY 


The in-plane x is isotropic and can be calculated as*”** 


1 2 
K = Koa = y È Ciba 
A 


(1) 


where V is the crystal volume, 4 is the phonon mode (comprising both 
a wave vector q and a phonon branch index j), C} is the heat capacity 
per mode, vj, is the group velocity of mode À along a direction, and 
Tja is the relaxation time. We use the nominal layer thicknesses h = 
4.20 A for both LB and DB silicene, corresponding to the van der 
Waals radii of silicon.” * Using the ShengBTE code,*°~>? we calculate 
the lattice thermal conductivity iteratively with the harmonic and 
anharmonic (third-order) interatomic force constants (IFCs) as the 
only inputs, which are obtained from first-principles calculations. 

We perform the calculations based on density functional theory 
(DFT) using the Vienna ab initio simulation package (VASP). We 
choose the generalized gradient approximation (GGA) in the 
Perdew—Burke—Ernzerhof (PBE) parametrization for the exchange- 
correlation functional. A plane-wave basis set with kinetic energy cutoff 
of 600 eV is employed. During structural relaxation for the unit cell of 
LB (DB) silicene, a 21 X 21 X 1 (15 X 15 X 1) q-mesh is used, until 
the energy differences are converged within 10~° eV, with a Hellman- 
Feynman force convergence threshold of 10 eV/A. A vacuum 
spacing larger than 15 A (20 A) is used for LB (DB) silicene to 
eliminate interactions between adjacent layers. 

By calculating the dynamical matrix through the linear response of 
electron density, the harmonic IFCs can be obtained based on density 
functional perturbation theory (DFPT).*° Using the supercell 
approach, a § X 5 X 1 (3 X 3 X 1) supercell with 5 x 5 x 1 (3 x 
3 X 1) q-mesh is used for LB (DB) silicene. The phonon dispersion 
relation is obtained using the harmonic IFCs, and subsequently the 
group velocity v,, and specific heat C} can be obtained. The third- 
order anharmonic IFCs play a significant role in determining the three- 
phonon scattering rate, which is the inverse of t,,. Using a supercell- 
based, finite-difference method, the anharmonic IFCs are calculated,*” 
and the same § X 5 X 1 (3 X 3 X 1) supercell is used for LB (DB) 
silicene with § x § X 1 (3 X 3 X 1) q-mesh. We include the 


interactions with the eighth nearest-neighbor atoms for both 
structures, which is well converged. 

The convergence of x with respect to q points is tested in our 
calculation. We introduce discretizationa of the Brillouin zone (BZ) 
into a [-centered regular grid of 90 x 90 x 1 (50 x 50 x 1) q points 
for LB (DB) silicene. A Gaussian function is used to enforce the 
conservation of energy in the three-phonon processes with a scale 
parameter of 0.5 for broadening.” 


3. RESULTS AND DISCUSSION 


3.1. Crystal Structures and Phonon Dispersions. Figure 
1 shows the optimized structure of LB and DB silicene. The 
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Figure 1. Top view and side view of (a) LB and (b) DB silicene. 


optimized LB silicene has a hexagonal structure with space 
group P3m1 and 2 Si atoms in one unit cell. A low-buckled 
configuration is found in Figure la, which is similar to LB 
germanene and stanene.”!°!*!!7 As shown in Figure 1b, the 
optimized DB silicene has a hexagonal structure with space 
group P6/mmm and 10 Si atoms in one unit cell. Two DB units 
are formed at the para positions of a honeycomb ring and the 
other Si atoms stay in the same plane, which is similar to DB 
stanene and DB Sn,Ge,.””*” The optimized geometries of LB 
and DB silicene are listed in Table 1. The lattice constant a and 


Table 1. Lattice Constant a, Buckling Height h, and 
Cohesive Energy E, of LB and DB Silicene® 


a (A) h (A) E. (eV/atom) 
LB silicene 3.87 0.44 4.77 
(3.87*°) (0.45*°) (4.57*°) 
(3.87*7) (0.44%) 
(3.88'*) (0.44'*) 
DB silicene 7.43 271 4.98 
(7.43'°) (2.71'°) 


“Other theoretical data are listed in parentheses for comparison. 


buckling height h of both LB and DB silicene are in good 
agreement with previous studies.'*'°*°*” In fact, when Si 
adatoms are placed on silicene, DB structures can be formed 
without any energy barrier,” and these adatoms gradually cover 
the surface uniformly.”* The DB structures can construct stable 
periodic structures with higher cohesive energy, because DB- 
DB coupling is attractive until the first nearest-neighbor 
separation,” 6 and the energy barrier to the migration of 
DBs along the straight path between the two first (second) 
nearest-neighbor atoms of silicene is 0.75 (0.95) eV.” As 
shown in Table 1, the calculated cohesive energy per Si atom 
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for DB silicene is 0.21 eV higher than that of LB silicene, 
indicating that DB silicene is more stable. 

High cohesive energy may not ensure that these two 
materials correspond to deep minima in the Born— 
Oppenheimer surface. Here, we investigate the dynamic 
stability of these two structures by the calculations of phonon 
dispersion in Figure 2. The calculated phonon dispersion of LB 
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Figure 2. Phonon dispersion for (a) LB and (b) DB silicene along 
different symmetry lines. 


silicene is in good agreement with previous works.*°-? The 
frequencies of all modes are positive over the whole BZ, 
therefore both LB and DB silicene are dynamically stable. 

As shown in Figure 2a, the longitudinal acoustic (LA) and 
transverse acoustic (TA) branches of LB silicene are linear near 
the I point, and out-of-plane acoustic (ZA) branch shows a 
quadratic trend with a linear component; while the ZA branch 
of DB silicene is totally quadratic near the I’ point, as shown in 
Figure 2b. This difference comes from the structural difference. 
DB silicene belongs to Dg, point group, whereas the buckled 
structure belongs to the point group symmetry of Dza which 
breaks the reflection symmetry, and therefore the Z modes of 
LB silicene couple with the XY modes.***” This results in the 
linear dispersion component of ZA phonons near the I’ point. 
Similar linear ZA dispersion has been observed in other buckled 
structures such as germanene and stanene.'”*”"*? 

The different shapes of ZA phonon dispersion lead to 
different group velocities v, for LB and DB silicene near the 
zone center. As shown in Figure 3a, a small but nonzero sound 
velocity for the ZA branch of silicene is observed in LB silicene, 
which is consistent with other studies.***° This is because the 
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Figure 3. Calculated phonon group velocities v, of (a) LB and (b) DB 
silicene in the whole BZ. 


dispersion of the ZA branch near the I point contains a linear 
component. The phonon group velocities of LB and DB 
silicene for the ZA, TA, and LA modes in the long-wavelength 
are listed in Table 2. It should be noticed that the sound 
velocities of LB silicene for the ZA, TA, and LA modes are 
larger than those of DB silicene. 


Table 2. Sound Velocities and Debye Temperature Op of LB 
and DB Silicene in Comparison with Graphene 


v,(ZA) (m/s) v,(TA) (m/s) »,(LA) (m/s)  @p (K) 


LB silicene 860 $570 9520 798.1 
DB silicene 0 4850 8700 753.4 
graphene 0 13 880°" 22 000°" 2300°° 


The Debye temperature ©p can be calculated with the 
maximum optical phonon frequency (Debye frequency) @, 


Oy = hwn / kpg (2) 


where h is the Planck constant and kg is the Boltzmann 
constant. The calculated ©p values of LB and DB silicene are 
given in Table 2. Because the Debye temperature reflects the 
magnitude of sound velocity,” it is expected that the larger the 
sound velocities, the higher the Debye temperature. The 
Debye temperature of LB silicene is higher than that of DB 
silicene, which is in consistency with the fact that the sound 
velocities of LB silicene are larger. For comparison, we also list 
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the phonon group velocities in the long-wavelength limit and 
Debye temperature of graphene in Table 2. The Op of 
graphene is much higher than that of LB and DB silicene. The 
Debye temperature is a measure of the temperature above 
which all modes begin to be excited and below which modes 
begin to be frozen out.°° Generally, low ©p corresponds to low 
thermal conductivity according to the aboved-mentioned 
Slack’s theory, >>" because lower Debye temperature indicates 
smaller phonon velocities and lower acoustic-phonon frequen- 
cies, and low acoustic-phonon frequencies further increase 
phonon populations, which subsequently increase phonon 
scattering rates.°°” 

The mode Griineisen parameters provide information about 
the anharmonic interactions and can be calculated directly from 
the anharmonic IFCs*” 


j (q*) oi 
= 1 apy Can Jeb, (a) iqR; 
Ifi = ~ 6a? X y X Z o o Tyry 
EET My 
(3) 


where a, P, and y are the Cartesian components, DF ay yr yr is 
the third-order anharmonic IFCs, d „(q*) and (q) are the 
phonon eigenvectors, and r,” p, is the position of n” atom in I” 
primitive cell. Figure 4 shows the mode Grüneisen parameters 
within the whole BZ for LB and DB silicene. The Grüneisen 
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Figure 4. Calculated Griineisen parameters for (a) LB and (b) DB 
silicene within the whole BZ. 


parameters of DB silicene are larger than those of LB silicene at 
low frequencies, indicating higher anharmonicity of DB silicene. 
Therefore, DB silicene has much stronger anharmonic phonon 
scattering than LB silicene at low frequencies. These large 
mode Griineisen parameters at low frequencies are crucial to 
the understanding of the phonon transport in DB silicene, 
because most of the heat is carried by low-frequency phonons. 

3.2. Thermal Conductivity. The room-temperature lattice 
thermal conductivity of naturally occurring (x,,,) and isotopi- 
cally pure (Kure) silicene are listed in Table 3. The calculated 


Table 3. Kaav Kpures P3, and Contribution of ZA, TA, LA, and 
Optical Phonon Modes toward the k,,, in All Studied 2D 
Group IV Crystals at 300 K 


Kyat K, ure 
(w (W ZA TA LA optical 
mK) mK) P; (%) (%) (%) (%) 
LB 27.72 28.85 0.15 38.98 21.63 20.97 18.42 
silicene 
DB 2.86 2.86 4.08 20.22 34.43 22.01 23.34 
silicene 


thermal conductivity of free-standing LB silicene is in good 
agreement with other theoretical results.**°*~°° Compared to 
LB silicene, much lower thermal conductivity is observed in DB 
silicene. It should be noticed that only free-standing silicene is 
considered here. However, in device applications, both LB and 
DB silicene are usually synthesized on various substrates, and as 
a result, the thermal conductivity will be further affected. In 
fact, for LB silicene, the influence of substrate on the thermal 
transport is in debate. Wang et al. have found that, in LB 
silicene, though the coupling of ZA modes to the substrate 
linearizes the phonon dispersion and increases the group 
velocity, the relaxation times of the ZA phonons are 
dramatically decreased, because the presence of the substrate 
breaks the reflection symmetry and introduces more phonon 
scattering." Thus, the presence of a SiO, substrate results in a 
large reduction to the thermal conductivity at 300 K. However, 
from nonequilibrium molecular dynamics simulations, Zhang et 
al. have found that by changing the surface crystal plane of the 
substrate, the thermal conductivity of LB silicene can be either 
suppressed or enhanced.*” 

The reduced thermal conductivity in DB silicene can be 
explained by conventional Slack’s theory. Lower Debye 
temperature in DB silicene results in decreased phonon 
velocities and higher phonon scattering rates since more 
phonon modes are active at a given temperature. °° The 
crystal structure of DB silicene is more complex: there are 10 
atoms in the unit cell of DB silicene, corresponding to 3 
acoustic and 27 optical phonon branches, indicating that a 
phonon can be scattered more easily in DB silicene (more 
three-phonon groups that can satisfy the energy and quasi- 
momentum conservations). In addition, large mode Griineisen 
paramters of DB silicene imply strong anharmonic phonon 
scattering as well. Consequently, the thermal conductivity of 
DB silicene is significantly reduced. Our results indicate that the 
lattice thermal conductivity of 2D group IV materials can be 
tuned by dumbbell structure formed from the addition of 
adatoms. 

As present in Table 3, the thermal conductivity of LB and 
DB silicene are insensitive to the naturally occurring isotope, 
since the Umklapp processes become dominant over isotopic 
scattering at high temperatures and consequently drive the 


thermal conductivity.°* In addition, DB silicene has much 
stronger anharmonic phonon scattering as mentioned above. 
Therefore, the isotopic influence on thermal conductivity is 
much less obvious in DB silicene than in LB silicene. 

Figure 5 presents the lattice thermal conductivity for LB and 
DB silicene with naturally occurring isotope concentrations as a 
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Figure 5. Thermal conductivity as a function of temperature for (a) 
LB and (b) DB silicene. 


function of temperature. The thermal conductivity decreases 
with increasing temperature, as expected for a phonon- 
dominated crystalline materials. We also compare the thermal 
conductivity calculated from the iterative and single-mode 
relaxation time approximation (SMRTA) solution. The differ- 
ence between the two approaches for DB silicene is much 
smaller than that for LB silicene. The SMRTA has no memory 
of the initial phonon distribution, because it assumes that 
phonons are excited individually. As a result, when the 
momentum-conserving Normal processes play an important 
role, the SMRTA typically does not work well.°° Our results 
indicate the dominant role of the Umklapp processes in 
phonon-phonon interactions of DB silicene over a wide 
temperature range. 

High TE performance requires the system to be a bad 
conductor for phonons. The ultralow lattice thermal con- 
ductivity of DB silicene (2.86W/mK at 300 K) is similar in 
magnitude to that of state-of-the-art TE materials such as PbTe 
(2.8 W/mK°°) and Cu,GeSe, (2.0—2.4 W/mK"). DB silicene 


may realize much higher thermoelectric efficiency considering 
the ultralow lattice thermal conductivity (2.86W/mK), which 
will induce further investigations on the performance of TE 
devices based on silicene. Compared to other 2D hexagonal 
materials, the lattice thermal conductivity of DB silicene is 
much lower than that of LB germanene (20.2 W/mK)! and 
stanene (11.6 to 26.2 W/mK),'”°' monolayer MoS, (30—100 
W/mK),°’—® and blue phosphorene (78 W/ mK), whereas 
it is at least 3 orders of magnitude smaller than that of graphene 
(3000—5000 W/mK).'***7! This is due to strong three- 
phonon scattering in DB silicene as we mentioned above. 

To understand the underlying mechanism of low lattice 
thermal conductivity of DB silicene, the key physical insight can 
be obtained by calculating the allowed phase space for three- 
phonon processes P}. In the temperature range where three- 
phonon processes are dominant, the total phase space for three- 
phonon processes P; is defined by?” 


} te 
B = Ze» + P ) 


30 (4) 
where Q is a normalization factor, > and 
p=) f daD (q) 
j (5) 
and 
Da =), fi dq’5(@(q) + w;(q’) 
ii" 
— w(q + q' - G)) (6) 


where Dj(q) corresponds to absorption processes, i.e., 
w(q)+@;(q') = @p(q+q'—G), whereas D7(q) corresponds 
to emission processes, ie, @(q) = @;(q’)+@j(q—q'—G). 
According to eqs 4—6, P, provides a measure of the phase space 
available for three-phonon processes of each phonon mode. 
Larger phase space for three-phonon processes usually implies a 
larger number of available scattering channels. Consequently, 
there is an inverse relationship between P, and the intrinsic 
lattice thermal conductivity of a material.°°”’* As shown in 
Table 3, DB silicene has larger P; than LB silicene, indicating a 
larger number of available scattering channels for phonons. 
This is consistent with the above analysis showing that the 
anharmonic phonon scattering in DB silicene is stronger. 

Furthermore, the relaxation time of each phonon mode for 
LB and DB ssilicene is extracted to reveal the scattering 
mechanism in detail. Figure 6 presents the isotopic relaxation 
time and the intrinsic three-phonon relaxation time of each 
phonon mode for these two materials. Isotopic scattering 
processes are frequency-dependent. The low-frequency pho- 
nons can transport nearly all the heat with little isotopic 
scattering, * thus at low frequencies we observe the relative 
long isotopic relaxation time of acoustic phonons. The total 
scattering rate is a sum of contributions from isotopic scattering 
and anharmonic three-phonon scattering, 6? 


1/7(q) = 1/7(q)"° + 1/7(q)"™ (7) 


The anharmonic three-phonon scattering is much stronger than 
isotopic scattering (corresponding to shorter relaxation time). 
Thus, with much shorter relaxation time, the three-phonon 
scattering dominates the total relaxation time. 

As shown in Figure 6a, the three-phonon relaxation time of 
LB silicene agrees well with previous results.” It should be 
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Figure 6. Isotopic and three-phonon relaxation time of each phonon 
mode for (a) LB and (b) DB silicene. 


noticed that 7;(q)*"" of LB silicene is much longer than that of 


DB silicene, leading to much higher lattice thermal con- 
ductivity. For DB silicene in Figure 6b, the enhancement of 
acoustic phonon scattering (corresponding to shorter relaxation 
time) is the main factor for the reduced lattice thermal 
conductivity. Moreover, the relaxation times of TA phonons in 
DB silicene are longer than ZA and LA phonons in the 
frequency range of 10—80 cm™. Therefore, the contribution of 
TA phonons to K,a of DB silicene is largest (34.43%), as shown 
in Table 3. 

3.3. Size Dependence of x. With ultralow lattice thermal 
conductivity, DB silicene may realize higher TE efficiency. This 
will motivate a systematic examination of the size dependence 
of thermal conductivity in these two materials for TE 
applications. To analyze the size dependence of thermal 
conductivity, we calculate the cumulative thermal conductivity 
with mean free path (MFP) below L by the following 
expression 


L<L 


K(L) = Kaa(L) = = È Cinta 

Vv" (8) 
The cumulative thermal conductivity at 300 K with respect to L 
for LB and DB silicene are shown in Figure 7a, b, respectively. 
For LB silicene, the accumulation function closely resembles a 
step function, with the major contribution by phonons with 
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Figure 7. Cumulative thermal conductivity for (a) LB and (b) DB 
silicene. 


MFP between 5 and 100 nm. For DB silicene, the major 
contribution comes from phonons with MFP between 0.5 and 
100 nm. 

Furthermore, when the size of sample gets smaller, the 
boundary scattering becomes dominant over other scattering 
mechanisms, and the x becomes proportional to a constant 
value lg which limits the MFP uniformly.” The Is can be 
calculated as 


lsa = K/Kgg (9) 


where Ksg is the thermal conductivity per Ig, in the small-grain 
limit. The Isg at 300 K is found to be 16.76 nm (2.34 nm) for 
LB (DB) silicene. These quantities are crucial to the 
manipulation of phonon transport in the small-grain limit 
such as nanowires. 

In both LB and DB silicene nanowires, the x will be further 
reduced with decreasing width due to stronger boundary 
scattering, as shown in Figure 8. At width about 13 and 4 nm 
for LB and DB ssilicene respectively, the lattice thermal 
conductivity drops about 50%, which means that x of DB 
silicene is less sensitive to nanostructuring size effects than LB 
silicene. Actually, higher lattice thermal conductivity arises 
typically because of weaker intrinsic phonon scattering, which 
causes phonons to have larger MFPs.” Consequently, these 
long-MFP phonons can be more easily blocked by boundaries 
of a similar size. Thus, because DB silicene has lower intrinsic 
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Figure 8. Thermal conductivities of (a) LB and (b) DB silicene 
nanowires along the [100] direction as a function of width. 


thermal conductivity, it is hard to reduce x further in 
nanostructures. 


4. CONCLUSION 


In summary, the lattice thermal conductivity of LB and DB 
silicene are calculated using first-principles calculations and the 
BTE for phonons. The x of both LB and DB silicene are found 
to be insensitive to naturally occurring isotopes. Much stronger 
anharmonic phonon scattering is observed in DB silicene as 
compared with LB silicene. Therefore, the predicted lattice 
thermal conductivity of DB silicene is much lower than other 
2D materials, which will induce further theoretical and 
experimental investigations on the TE performance for these 
two structures of silicene. Finally, for designing nanostructures 
in TE applications, we investigate the size-dependent thermal 
conductivity in detail. 
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